Graviton Mass from Close White Dwarf Binaries Detectable with LISA 
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The arrival times of gravitational waves and optical light from orbiting binaries provide a mechanism 
to understand the propagation speed of gravity when compared to that of light or electromagnetic 
radiation. This is achieved with a measurement of any offset between optically derived orbital 
phase related to that derived from gravitational wave data, at a specified location of one binary 
component with respect to the other. Using a sample of close white dwarf binaries (CWDBs) 
detectable with the Laser Interferometer Space Antenna (LISA) and optical light curve data related 
to binary eclipses from meter-class telescopes for the same sample, we determine the accuracy to 
which orbital phase differences can be extracted. We consider an application of these measurements 
involving a variation to the speed of gravity, when compared to the speed of light, due to a massive 
graviton. For a subsample of ~ 400 CWDBs with high signal-to-noise gravitational wave and optical 
data with magnitudes brighter than 25, the combined upper limit on the graviton mass is at the level 
of ~ 6 X 10"^'' eV. This limit is two orders of magnitude better than the present limit derived by 
Yukawa-correction arguments related to the Newtonian potential and applied to the Solar-system. 



I. INTRODUCTION 

In a recent paper, we discussed the optical foUowup 
study of close white dwarf binaries (CWDBs) that will 
be detectable with the Laser Interferometer Space An- 
tenna (LISA) mission [1]. Recent estimates suggest that 
close to 3000 binaries will be detected via gravitational 
waves (GWs) at frequencies above 3 mHz, with most of 
them restricted to the frequency range between 3 mHz 
and 6 mHz [2-4] . These binaries will only be localized to 
the accuracy of ~ one sqr. degree with a three-year LISA 
observation based on the Doppler effect related to the or- 
bit around the Sun and the amplitude modulation due to 
the rotation of LISA detectors. For a precise location of 
an individual binary, optical observations are required; 
an easy identification of the binary is facilitated by the 
presence of eclipses in the optical light curve with a pe- 
riod twice the GW period. Considering the impact on 
astronomy and astrophysics, in general, one can certainly 
expect detailed optical followup observations of LISA bi- 
naries [1]. In addition to CWDB related physics, such 
as tidal heating, using both the optical light curve and 
gravitational wave data, we can also extract important 
information related to certain fundamental parameters 
in physics. We discuss such a possibility in the present 
paper by considering observational prospects in the LISA 
era. 

In addition to a localization of the direction from which 
GWs are emitted, with gravitational waves, LISA also 
allow a measurement of the binary distance, the chirp 
mass, the binary orbital inclination angle, and the binary 
orbital phase (modulus tt). The optical followup obser- 
vations also allow a measurement of the orbital phase. 
In general, one expects the orbital phase measured with 
gravity wave data to agree precisely with that of the opti- 
cal light curve. A difference in the two is only expected if 
optical light travel at a different speed than that of grav- 
ity, which determines the propagation of gravitational 



waves. While a difference is not expected, a potential 
non-zero mass related to the graviton particle can lead 
to a slight difference. Thus, a comparison of the or- 
bital phase related to an optical eclipsing light curve of 
CWDBs with the phase determined from LISA data will 
allow a constraint on the graviton mass. 

Note that previous studies have considered such a pos- 
sibility, again with respect to LIGO and LISA data [5-8] . 
While one of the considered methods involve the disper- 
sion of the gravity waves alone, as a function of the fre- 
quency [5,6], the comparison between optical or electro- 
magnetic data and gravity waves is expected to provide 
an improved constraint. Previous analyses, however, only 
considered a handful of objects which are already known 
to be eclipsing binaries in optical data and are expected 
to be GW sources detectable with LISA (e.g., 4U1820- 
30) [8]. The present study considers a sample of CWDBs 
detectable and locatable with LISA above its confusion 
noise and is bright enough optically for followup studies 
with a few meter-class telescope. 

The paper is organized as follows. In § II, we briefly 
discuss observations with LISA, requirements for an opti- 
cal followup study and, then, consider the extent to which 
orbital phases of individual binaries can be determined. 
In § HI, we put these measurements in the context of an 
improved limit on the graviton mass. We conclude with 
a summary in § IV. 



II. CLOSE WHITE DWARF BINARIES 

First, we will review CWDB detections with LISA in 
gravitational waves and then move on to discuss aspects 
related to the optical light curve. Our focus here would 
be related to the phase measurement related to the bi- 
nary orbit, while issues related to the CWDB localiza- 
tion with LISA data and optical followup observations 
are discussed in Rcf. [1], following initial calculations in 
Refs. [9,13,14]. 
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A. LISA Observations 

Briefly, with gravitational waves, one observes two 

components given by the quadnipolc approximation in 
the principle polarization coordinate [10] 



= A cos 
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where n, given by {9s, (ps), and L are unit directional 
vectors to the binary, from the observer, and the angu- 
lar momentum of the binary, respectively. Here, ip{t) is 
the phase resulting from the Doppler phase due to the 
revolution of LISA around the Sun: 



ip{t) = 27r/i?sin^5 cos 
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where R = 1 AU and (po are integral constants. The 
principle polarization coordinate is determined by two 
orthonormal vectors X and Y that are normal to the 
direction of the source n through 
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CWDBs are expected to have circular orbits due to 
the tidal interaction in their early evolutional stage. The 
amplitude of the wave is given by 
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where D is the distance to the GW source. Note that 
the GW frequency (= 2/Porb where Porb is the orbital 
period) for a circular orbit is related to the total mass 
Ml + M2 and the separation a of the binary via 
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while the time variation of this frequency is 
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With gravitational waves, one can estimate a total of 
8 independent quantities: A, /, /, (po, location (n, 2 pa- 
rameters), and the direction of the angular momentum, L 
(2 parameters). At higher GW frequencies, one can also 
extract meaningful information on the second derivative 
of the frequency /. The orbital inclination angle is given 
by cos~^(n • L). In terms of physical quantities of inter- 
est, with A, /, /, one extracts D and M; in addition to 
these parameters, one can also constrain a combination 
of Mtot and a with Eq. 5. Note that relations (4) and 
(6) are given for Newtonian point particle systems as the 
correction related to the finite size of the binary is not 
significant. 

With parameters related to h and L, one can recon- 
struct the orbital phase, say at a certain time t, which 
we define as 
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when the chirp mass is given by M = Mf^^ M2^^ / {Mi 
M2)^/^ and the total mass Mtot = Mi + M2. 



(p{t) = tan"^ 
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where Z = Lx X and f{t) is the unit vector between the 
two masses. This vector changes as a function of time as 
the binary components orbit with respect to a common 
center of mass. 

For the present discussion, we extend the calculation 
in Ref. [1] to study how well the orbital phase in each of 
LISA detectable CWDBs can be determined. In order to 
generate the necessary representation of LISA observa- 
tions, we make use of the sample studied in [1] based on 
population synthesis code used by Ref. [11]. This binary 
sample is distributed in the galaxy following 



p{R, z) = posech 
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where Zq = 200 pc and Rq = 2.5 kpc [12]. For the 
sample of ~ 3000 binaries, we now randomly assign an 
orientation L. This means that the inclination cos i is 
distributed uniformly in < cosi < 1. For an optical 
eclipse to be observable, we require that the cosine of 
the inclination angle be less than a minimum such that 
cosz < cosiniin = (-Ri + R2) / <i- For simplicity, here, we 
assume that binaries are equal mass. For a given or- 
bital frequency and a chirp mass, we calculate the ratio 
(i?i + R2)/a as a function oi q = Mi /Mid using a typi- 
cal mass-radius relation of white dwarfs given in Ref. [15] . 
We found that the ratio -f i?2)/a takes a minimum 
value at q = 0.5 and increases only ~ 10% at q = 0.3 
or q = 0.7. This fractional change depends very weakly 
on the chirp mass. Considering the facts that the mass 
ratio M1/M2 for a CWDB is expected to be close to 1 
[16], our equal-mass assumption is conservative and safe. 
^5pnder these parameters, we extracted a siibsamplc of 
objects from the total original sample that is expected to 
emit detectable GWs with LISA and also show eclipses 
in the optical light curve. Since a typically is few times 
+R2), the probability to observe an eclipse is rather 
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high and this subsamplc amounts to about 1400 objects 
with an overaU probabihty of ^ 0.45. The subsample size 
is likely to vary higher by up to ~ 100 objects when our 
assumption related to equal mass binaries is dropped. 

The distribution of phase errors for this sample of ~ 
3000 CWDBs detectable with LISA is shown in Fig. 1. 
The phase error is calculated using the Fisher matrix 
approach discussed in Ref. [13], where all 8 parameters 
related to these observations arc derived simultaneously 
from GW data. As in Ref. [13], we use expressions related 
to LISA noise appropriate for the long wavelength limit, 
when compared to arm-length, with a transfer function 
for the finite size of the arm-length that subsequently 
corrects for the first assumption. This prescription is 
very effective for nearly monochromatic sources [17]. The 
noise curve in Ref. [18] is used here. We consider three 
cases with a total observational duration with LISA of 
1, 3 and 5 years. The lower set of lines show the phase 
error related to the subsample that is expected to show 
eclipses in the optical light curve, simply selected based 
on the inclination angle. This subsample can be optically 
identified and localized precisely; however, assuming that 
the location is exactly known is not important when the 
observational duration is more than 2 years (see, e.g., 
Ref. [13] for details on parameter improvement when the 
binary location is assumed to be known precisely a pri- 
ori). 

As discussed in Ref. [8] , the phase error related to GWs 
is simply given by 

'5'/'GW = f (I) \l + 0{S/N)-'] (9) 

where, a is order unity or more and accounts for the fact 
that more than one parameter is to be determined from 
GW data. The factor of 2 in above accounts for the fact 
that we are determining the phase error related to the 
binary orbit, which is 1/2 the gravity- wave phase. The 
phase errors are generally at the level of ~ 0.1, when 
the observational duration is more than 3 years. This is 
substantially smaller than what was previously suggested 
in the literature (c.f. [7])* 

In Fig. 1, for reference, we also show the phase error 
related to the nearby CWDB sample (at distances below 
8 kpc). There are roughly 1000 CWDBs and we expect 



*The phase error is independent of the timing error and the 
extent to which the period is known, as assumed in estimates 
by [7] . This can be simply understood by Fourier decompos- 
ing the hght curve. The phase is simply related to the angle 
between that of a Fourier mode and an arbitrary vector, re- 
lated to the location from which the phase is measured, and 
the extent to which this angle can be determined depends on 
the length of the Fourier mode, which in return is determined 
by the signal-to-noise of the light curve and not the timing 
accuracy. 



roughly ^ 500 of them to show eclipses that can be de- 
tected in the optical light curve. The phase error related 
to the nearby sample is better than the whole sample 
and results from the fact that the signal-to-noise for an 
individual detection is higher given the close distance. 

Note that the orbital phase determination is substan- 
tially improved with observational durations, Tobs, at the 
level of 3 years. While one generally expects the phase 
error determination to scale as vTobs, there is a substan- 
tial improvement between one and three-year durations 
beyond that implied by this simple scaling [13] . This in- 
crease resiilts from the fact that significant degeneracies, 
between the binary location and other parameters, are 
broken with observational durations greater than 2 years. 
We illustrate this in Fig. 2 where we show the phase error 
distributions for parameter determination that involve (i) 
all eight parameters, (ii) six parameters with location as- 
sumed to be known, and (iii) four parameters with both 
location and orbital orientation assumed to be known. 
The 1 year curves with cases (i) and (ii) show the extent 
to which degeneracies affect the orbital phase determina- 
tion. With location assumed to be known, for an observa- 
tional duration of 1 year, the degeneracy is broken and a 
substantial improvement in the orbital phase determina- 
tion is achieved. In contrast, curves for cases (i) and (ii) 
almost overlap when Tobs =3yrs or 5yrs. These results 
are instructive to understand how the impact of precise 
localization (through optical identification) on parame- 
ter estimation depends with the observational duration 
of LISA. 

Considering the observational situation relevant for 
this paper, we can implicitly assume that individual bi- 
nary locations are determined precisely by the identifica- 
tion of optical counterparts. However, when we consider 
observational durations of order 3 years or more for GWs, 
the final graviton mass limit we calculate in the present 
paper is essentially independent of whether the binary 
location is assumed to be known or not in the GW phase 
measurement. One can obtain an improvement of a fac- 
tor of 2 in the GW orbital phase when the orientation of 
the binary is also assumed to be known a priori. We do 
not pursue this possibility, however, since this involves 
detailed modeling of the optical light curve or obtaining 
additional data at different wavelengths of the electro- 
magnetic radiation spectrum. For example, we might get 
limited information on the inclination angle by modeling 
the eclipse shape, but the polarization angle would not 
be determined. 



B. Optical data 

The optical magnitude distribution of the CWDB bi- 
nary sample that is expected to be localized with LISA 
is discussed in Ref. [1]. Most optical magnitudes are at 
the level of 25 and fainter. Thus, we will only restrict the 
discussion to the subsample that is bright. Similarly, the 
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extent to which the phase angle can be determined from 
the optical light curve depends on the single-to-noise ra- 
tio of the optical light curve, and we write 

5<^opt=/3(|) [1 + 0{S/N)-'] , (10) 

where now /3 is again a parameter of order unity or more 
as we are trying to determine additional parameters from 
the optical data as well. Under the assumption that the 
period is a priori known, from the GW data, and the as- 
sumption that the time derivative of the period is zero, 
which is a very safe assumption since P ~ 10~ ^"ss"^ 
for binaries at 3 mHz with chirp mass ~ 0.45 Mq, /3 is 
generally at the level of 1 or slightly above. Since the 
observational duration of GWs, Tobs^GW ~ 3 years or 
more, is larger than the duration of optical observations 
(~ 100 minutes), our assumption of neglecting P for op- 
tical analysis if safe, though, we include P, through /, in 
estimating parameters related to GWs. 

For a meter-class telescope, with an integration time 
of tint, the expected signal-to-noise ratio can be written 
as 

\/ "SwD^int + iSskyiint + "Sdet^int 

where S'wd, S'sky and S'dot arc the number of source, sky 
background and detector photons, respectively, per unit- 
time. For a star of magnitude m, imaged with a telescope 
of diameter D, efficiency e, the number of photons, in a 
second, is 




(12) 



In the limit that sky and detector backgrounds are not 
important, the signal-to-noise, for a 25-second integra- 
tion is 




(13) 

In reality, at these faint magnitudes, the signal-to-noise 
ratio is, however, not determined only by the flux from 
the object, but rather by the sky and the detector back- 
ground, and the true signal-to-noisc, for a 25 second ob- 
servation, is at the level of 76 for m = 20, 12 for m = 22.5 
and ~ 1.5 for m = 25, in the case of the 4-m Cerro-Tololo 
and Kitt Peak telescopes of the National Optical Astro- 
nomical Observatories'^. In our calculations, we use these 



numbers, but restrict our attention to the subsample at 
a distance below 8 kpc and with magnitudes brighter 
than 25; with this additional restriction, the subsam- 
ple reduces to about ~ 400 CWDBs, for which reliable 
phase measurements may be possible from both GWs, 
with LISA, and optical data. 

Given the magnitude distribution of the CWDBs that 
are expected to show optical eclipses, we can estimate the 
expected distribution of errors in the orbital phase, de- 
termined from electromagnetic radiation. We summarize 
our results in Fig. 3 for the subsample of CWDBs, where 
we also include extinction [19], following the modeling in 
Ref. [1]. Most magnitudes arc at the level of 25 and be- 
low, though for the sample with distances below 8 kpc, 
magnitudes are at the level of 22, albeit a broad distribu- 
tion over the range of 20 to 25. While magnitudes here 
are simply calculated based on white dwarf cooling, given 
the age distribution determined in Ref. [1], there is an ad- 
ditional effect that can potentially brighten CWDBs due 
to tidal heating by one white dwarf on the other due to 
their close separation [20]. The magnitude distribution 
gains few magnitudes as shown in Fig. 3, but we do not 
consider this possibility given the complex physics associ- 
ated with tidal heating such as the possibility that tidal 
heating is not uniform and will only lead to hot-spots 
among others. However, we note that, from detailed fol- 
lowup observation of nearby CWDBs, one can obtain in- 
formation that would potentially affect the orbital phase 
measurement from optical data, e.g. brightness distribu- 
tion of the surface of white dwarf. Note also that inter- 
acting binaries, such as, AM CVn systems would be more 
luminous, but they generally have smaller chirp masses 
(weaker gravitational wave signals), and estimation of 
their orbital phases from optical data would not be as 
simple as CWDBs (see e.g. [7]). 

III. GRAVITON MASS 

With phase errors related to the GW and optical light 
curve data, one can consider a combined study involving 
the propagation speeds of gravity (cg) and light (c). The 
first obvious application is to determined if Cg = c and 
to study any departures from it. Compared to many 
other previous tests, this is an ideal scenario since one is 
comparing directly the gravitational radiation and with 
optical light with a calibrated zero-point with respect to 
a certain location of the orbit. The presence or the lack 
of any departure in the optically derived orbital phase 
from that of GWs can then be used to constrain some 
aspects related to the propagation that potentially make 
two speeds depart from each other. 

As such a possibility, we consider the presence of a 
massive graviton such that the speed of gravity is modi- 
fied 
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FIG. 1. The phase error distribution for GWs detectable 
with LISA for a sample ~ 3000 CWDBs, and assuming a to- 
tal observational duration for 1, 3 and 5 years. The error 
in orbital phase is calculated following the Fisher matrix ap- 
proach, with the total of eight parameters to be determine 
from gravity-wave data, including binary location. The lower 
curves are the error in orbital phase for the CWDBs that 
are expected to show eclipses. For comparison, we also show 
phase errors related to the subsample, about ~ 500 objects, 
that are at distances below 8 kpc and are expected to show 
echpses in the optical light curve, under the condition that 
cosi < {Ri + R2)/a, assuming a total observational dura- 
tion of 5 years. On average, the phase measurement with the 
nearby sample is better since the expected signal-to-noise for 
gravitational wave detection is higher. 



2 4- 



1/2 



1 - 1/2 



? 4 

rrigC 



(14) 



For the orbital binary, assuming a gravity wave frequency 
of /(= 2/Porb), and at a distance D the two phases are 



fopt 
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Simplifying with Eqs. 14 and defining A = 0gw ~ 0opt, 
one leads to the constraint that 



nigc^ = hf 



TTfD 
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(17) 



In the absence of a measurable phase difference A = 0, 
an upper limit on rrig can be obtained with the limit on 
A such that 
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FIG. 2. We show the role of prior knowledge related to 
the binary location and orientation of when the orbital phase 
errors is determined from GWs detectable with LISA for a 
sample ~ 3000 CWDBs, and, again, assuming a total ob- 
servational duration of 1 (solid-lines), 3 (dotted- lines) and 5 
(dashed-lines) years. The three sets of curves with triangles, 
squares and circles, represent the case with (i) no prior knowl- 
edge is available and all eight parameters are to be determined 
from the GW data, (ii) location assumed to be known, and 
(iii) both location and orientation of the binary is assumed to 
be known, respectively. In general, the location and other pa- 
rameters are highly correlated and precise localization leads 
to a substantial reduction in the accuracy of phase determi- 
nation for observational durations of two years and less. For 
observational durations greater than 3 years, the improve- 
ment in phase error related to a priori known location is not 
significant when compared to the case with all parameters, in- 
cluding location, are to be determined from the data. Curves 
for cases (i) and (ii) are almost overlapped for Tots =3yrs and 
5yrs. Additional degeneracies are broken when the binary ori- 
entation is also assumed to be known and results in a better 
determination of orbital phase. We do not pursue this possi- 
bility since it involves detailed modeling of the optical light 



In the limit that S(/)gw >> <^^opt, the limit on graviton 
mass can be written in terms of the signal-to-noise to 
which LISA detects gravitational waves [8] 
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and since the GW signal-to-noise simply scales as l/D, 
this limit is independent on the distance at which the 
binary is located [8] . In contrast the signal to noise ratio 
of electro-magnetic waves scales as l/D^, such that in 
the opposite limit, the limit on graviton mass is no longer 
independent of the distance. 

As shown in Fig. 1, and implied from Fig. 3, either 
one of these assumptions, S(/)gw >> <^0opt or S<pGW << 
(18) (5(/)opt, however, is not correct since the error to which 
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FIG. 3. The V-band magnitude distribution for the CWDB 
sample detectable with LISA and for which phase measure- 
ments may be possible. We have taken into account the spa- 
tial distribution in addition to extinction (solid line). The 
long-dashed line show the subsample at distances below 8 
kpc. This subsample, on average, is brighter given the nearby 
distance and relatively low extinction, though the distribution 
is some what broad. The dash-dotted line shows the possible 
effects of the tidal heating with using the extreme model in 
[20]. 



both optical and gravity-wave phases are determined is 
at the same order. Thus, one should include both errors 
in estimating the hmit on graviton mass. For the present 
discussion, we make use of the sample at a distance be- 
low 8 kpc. The subsample is brighter and is within the 
magnitude limits comfortably reachable with few meter 
or more-class telescopes that are currently available. We 
make use of a fraction of 400 objects with magnitudes 
brighter than 25, and with distances below 8 kpc. Here, 
we neglect the uncertainty in the distance estimation re- 
lated to LISA data since its error is in the second or- 
der. The individual distances are, however, statistically 
included by accounting for the distance distribution func- 
tion based on our model assumption related to the spatial 
distribution following Eq. 8. 

Note that the observations we consider are indepen- 
dent of the localization of CWDBs since this can be 
achieved at a high signal-to- noise ratio by imaging at time 
intervals corresponding to expected minima and maxima 
in the optical light curve based on prior phase informa- 
tion from the gravity-wave data. To obtain the optical 
phase precisely, detailed sampling of the light curve is re- 
quired and will involve imaging at time intervals of order 
25 seconds, given that the periods are mostly between 
400 to 800 seconds. The orbital phase errors related to 
optical and GW data are summarized in Fig. 4 for the 
subsample of CWDBs. The phase errors are generally at 
the level of 10^^ for GWs, while the orbital phase error 
from the optical data follow a broad distribution. The er- 
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FIG. 4. Top: The error in orbital phase distribution for a 
subsample of ~ 400 CWDBs at distances below 8 kpc and 
at magnitudes brighter than 25. We show separately the or- 
bital phase error from GW and optical data and the combined 
error related to the measurement of the phase difference in 
each of these cases. For each binary, we assume optical data 
of over 10 orbital cycles or roughly over a time interval of 
^ 6000seconds, and consider the signal-to-noise measures as 
expected from a four- meter class telescopes (see, equations 11 
to 13 and the discussion below Eq. 13). For the whole sample 
of 400 CWDBs, such a monitoring project takes a total of ~ 
90 days assuming 8 hours of useful observational time each 
day. Bottom: the distribution of graviton mass limits derived 
from the above sample, together with its distance and fre- 
quency distributions. The limit, for an individual CWDB, is 
generally at the level of 10~^^ eV, but for the sample as a 
whole, this limit improves to about ~ 6 x 10^^* eV. The dot- 
ted line is the results only with considering the phase errors of 
gravitational waves. The combined limit becomes ~ 3 x 10"^^ 
eV. 



ror related to the combined phase difference is also shown 
in Fig. 4, which is calculated using physical properties of 
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the subsamplc. In general, note that the phase error 
distribution related to GWs and optical data are not in- 
dependent; the CWDBs with better determined orbital 
phases with GW data are expected to be nearby and, 
thus, to be bright optically; For these nearby ones, or- 
bital phases from optical light curve data are also better 
determined. 

In Fig. 4 bottom plot, we summarize the extent to 
which the graviton mass can be constrained from each of 
the CWDBs at distances below 8 kpc and with magni- 
tudes better than 25. The one-sigma limit on the gravi- 
ton mass is generally at the level of 10~^^ eV. If phase 
information related to the whole sample can be combined 
for a single constraint as whole, the combined limit, by 
taking {J^i improves to about 6 x 10"^'* eV 

at the one-sigma level. This can be compared to cur- 
rent limits on the graviton mass. Based on the Yukawa- 
correction to a Newtonian potential, in the presence of a 
massive graviton, the limit from Solar-system dynamics 
is 4.4 X 10^^^ cV [21], though the limit improves sub- 
stantially to 2 X 10"^^ eV with galaxy clusters at Mpc 
size-scales [22,23]. The 90% confidence limit from the 
orbital decay of two- binary pulsars, again an indirect es- 
timate on the graviton mass, is 7.6 x 10^^" cV [24]. The 
direct limits from GW and optical data are discussed in 
Refs. [7] and [8] ; these limits are generally above the level 
we have discussed, as they usually involve either a sin- 
gle binary or different assumptions and estimates. Our 
limit, for the sample as a whole, is consistent with the 
best possible limit from LISA data implied in Ref. [8]. 

In general, we believe that the CWDB sample related 
to LISA provides a well controlled sample with adequate 
optical followup opportunities to obtain a reliable limit 
on the graviton mass as a whole. When obtained, such 
a limit will be independent of indirect assumptions and 
will be directly based on the propagation speeds of grav- 
ity and light. Finally we comment on a potential con- 
taminant related to this measurement involving gravi- 
tational lensing of the background binaries by the fore- 
ground mass distribution. Note that the lensing prob- 
ability in our galaxy has been measured to be at the 
level of T ~ 10^^ [25], suggesting that the probabil- 
ity to observe a lensed CWDB is substantially small. 
Since gravitational waves from CWDBs have much lower 
frequencies 10~^Hz) than optical waves ~ lO-'^^Hz, 
the critical frequency for lensing with a mass Miens is 
U - (GM,,„,/c3)-i ~ lO4(Mie„./M0)-iHz. For ex- 
ample, the amplification effect is almost negligible for 
waves with / ^ fcr where the geometrical optics ap- 
proximation is not valid [26-28]. There might be ef- 
fective time delay between the gravitational wave and 
the long wave limit (electro-magnetic waves) at most 
GMipns/c^ ^ 10~^seconds with a phase shift at the level 
of ~ 10~^ radians. This, however, is much smaller than 
our resolution A^/(27r/) ~ 10 seconds and well below 
the phase measurement error. Given the low probabil- 
ity for lensing and the small phase difference expected, 
we do not consider lensing contamination to be a major 



source of concern. 

Other potential contaminants include modifications to 
the speed of light, relative to gravity waves, based on 
refractive fluctuations through the neutral interstellar 
medium, such as gas clouds. Gravity waves will pass 
such clouds unobscured, while density gradients within 
the medium will lead to refraction and, thus, modifi- 
cations to the path length of the optical light. Here, 
again, expectations are relatively small modifications, if 
any, but we suggest that further studies may be necessary 
given the importance of combined GW and optical stud- 
ies in the LISA era and the unique sample of CWDBs 
that LISA data are expected to provide. 

IV. SUMMARY 

The arrival times of gravitational waves and optical 
light from orbiting binaries provide a mechanism to un- 
derstand the propagation effects of gravity, when com- 
pared to light. This is easily achieved via binary orbital 
phase measurement and by looking for an phase offset in 
the optically derived orbital phase related to that derived 
from gravity wave data. Using a sample of close white 
dwarf binaries (CWDBs) detectable with the Laser In- 
terferometer Space Antenna (LISA) and the associated 
optical light curve data related to binary eclipses, we de- 
termine the accuracy to which orbital phase difference 
can be extracted. We consider an application of these 
measurements, determining an upper limit on the gravi- 
ton mass. 

For a subsample of ~ 400 CWDBs with high signal- 
to-noise gravity wave and optical data with magnitudes 
brighter than 25, the combined upper limit on the gravi- 
ton mass is at the level of ~ 6 x 10~^^ eV, which is 
two orders of magnitude better than the limit derived by 
Yukawa-correction arguments related to the Newtonian 
potential and applied to the Solar-system. 
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